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We investigate an extended ±J Ising spin glass model by using a gauge symmetry. This 
model has ± Ji interactions and ± J 2 interactions. We show that a gauge symmetry is usable 
to study this model. The exact internal energy, the rigorous upper bound of the specific 
heat and some rigorous relations for correlation functions and order parameters are shown 
by using the gauge symmetry. The results are rigorous, and do not depend on any lattice 
shape. A part of our results, e.g., the value of the exact internal energy should be useful for 
checking the computer programs for investigating this model. In addition, we find that the 
present solutions are general solutions which include the solutions on the so-called Nishimori 
line for the conventional ±J Ising spin glass model and the solutions for the bond-diluted 
±J Ising spin glass model. 

The theoretical studies of spin glasses have been widely doneP^® For the 
universality classes of spin glasses at zero temperature on the square lattice, the 
±J Ising spin glass model indicates that the stiffness exponent 9 is zero for the 
symmetric distribution of randomness an d the Gaussian Ising spin glass model 
indicates 9 ~ —0.28 for the symmetric distribution of randomness.^'^ The stiffness 
exponent 9 is investigated from the system-size dependence of domain wall energy.^ 
The present model, i.e., an extended ±J Ising spin glass model has ±Ji interactions 
and ±J2 interactions. When J2/J1 = 2, the 9 is zero for the symmetric distribution 
of randomness.^ Therefore, when J2/J1 = 2, the model has the universality class 
of the conventional ±J Ising spin glass model. On the other hand, when J2/J1 is 
(V5 + l)/2 (the golden mean) ~ 1.618, the model is called the irrational model, and 
9 ~ —0.29 for the symmetric distribution of randomness is concluded^ Therefore, 
the universality class of the irrational model is associated with the universality class 
of the Gaussian Ising spin glass model rather than that of the conventional ± J Ising 
spin glass model. It is pointed out that this difference between the universality classes 
depends on the difference whether the energy level in the histogram for the domain 
wall energies is quantized or continuous!^ This property of the present model for 
universality class is remarkable. The studies for this model are seen in Refsj2]),|3]). 
We investigate this model by using a gauge symmetry. 

Analysis by using gauge symmetries makes rigorous arguments possible. How- 
ever, the applicable models are not so many. Here, the applicable models are models 
that physical quantities and/or relations are obtained by using gauge symmetries. 
As the applicable models, there are the conventional ± J Ising spin glass modelP® 
the bond-diluted ± J Ising spin glass modelP the Gaussian Ising spin glass model 
the Potts gauge glass model 8 ^ and the XY gauge glass modeP 1 '® for example. We 
show that a gauge symmetry is usable to study the present model. In the applica- 
tions of many other methods, a lattice shape is supposed in advance, and the results 
are calculated on the lattice. On the other hand, in this method using a gauge 
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symmetry, any lattice shape is not supposed in advance. 

The Hamiltonian for the Edwards- Anderson Ising spin glass model, T~L, is given 

% = — JijSiSj , (1) 

(id) 

where denotes nearest-neighbor pairs, Si is the state of the spin at a site i, and 
Si = ±1. \Jij\ is the strength of the exchange interaction between the spins at sites 
i and j. The value of Jjj is given by a distribution P(Jjj) which is^ 1 

P(Ji,j) =pS{Jij - J) +qS(J iij + J) + r 5(Jij - a J) +u5{J lJ + a J) , (2) 

where J > 0, and 5 is the Dirac delta function, p, q, r and u are probabilities. 
p is the probability that the interaction has the strength of J and the interaction 
is ferromagnetic, q is the probability that the interaction has the strength of J 
and the interaction is antiferromagnetic. r is the probability that, if a > 0, the 
interaction has the strength of a J and the interaction is ferromagnetic, r is the 
probability that, if a < 0, the interaction has the strength of \a\J and the interaction 
is antiferromagnetic. u is the probability that, if a > 0, the interaction has the 
strength of aJ and the interaction is antiferromagnetic. u is the probability that, if 
a < 0, the interaction has the strength of \a\J and the interaction is ferromagnetic, 
r + u is the probability that, if a = 0, the interaction is diluted. For p, q, r and u, 
we have 

p + q + r + u = l. (3) 

a is a real number. When a = 1 or —1, the model is the conventional ±J Ising spin 
glass model. When a = 0, the model is the bond-diluted ± J Ising spin glass model. 
When a = (\/5 + l)/2 or (\/5 — l)/2, the model is the irrational model mentioned 
in Ref.ED. 

To calculate thermodynamic quantities, a gauge transformation is used and is 
performed by"™™''"*^' 

OiOj , Si -4 S'jCJj , (4) 

where Ui is a variable at a site i, and takes 1 or —1 arbitrarily. For the conventional 
Ising spin glass model, this gauge transformation has no effect on thermodynamic 
quantities.^ For the present model, also this gauge transformation has no effect on 
thermodynamic quantities. By using the gauge transformation, the T~L part becomes 

By using Eq. ([2]), the distribution P(Jij) is rewritten as 

P(J itj ) = Ae^ ) <J + ^ , (5) 

A e /B ( p ) J 2 +l3pJ + e l3 ( jPj 2 -l3pJ + e /3^ 2) a 2 J 2 +/3 P aJ + e /3^ 2) a 2 J 2 -/3 P aJ ' ^ 

^e^ )j2+ ^ J =p, (7) 
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Ae ^ )j2 -^ J = q, 

Ae P { p ] a 2 J 2 +p P aJ =J . 
Ae i3^a 2 J 2 -f} P aJ = u _ 



(2) 

Pp , ftp and a are respectively 



i , ( pq 

2(1 -a 2 )J 2 U \ru 



ln(r/u) 



ln(p/g) ' 

By using the gauge transformation, the distribution P(Jij) part becomes 



(8) 
(9) 
(10) 



(11) 
(12) 
(13) 
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(14) 



where iV is the number of sites, and Np is the number of nearest-neighbor pairs in 
the whole system. 

By using Eqs. ([I]) - (|14p . when the value of f3p is consistent with the value of 
the inverse temperature /3, the exact internal energy, [{ / H)t p ]r, is obtained as 



W) 



T P \R 



e e n nji, m ) 

{Ji.m} (l,m) 
A N B 



,(2) ■ 



Pp E(!,m> Jl,niSlS„ 



2 N 

-JN B 



E E(- J «^'^) e ^ E<wJ! - ro+/3pE<! ' ra> ' 7i ' m5i5T ' 

<i,i)«,m}{^} 

ln(r/u) 



p — q + 



ln(p/g) 



(r — u) 



(15) 



where ()t p is the thermal average at the temperature T = Tp, []p is the random 
configuration average, ftp = 1/kpTp, and kp is the Boltzmann constant. The value 
of the exact internal energy should be useful for checking the computer programs for 
investigating this model for example. 
The specific heat C is given by 



C = k B p 2 (l(H 2 ) T ] R -[(H) T } R ) 



(16) 
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where ( )t is the thermal average at the temperature T. By performing a similar 
calculation with the calculation in Eq. (fl"5j) . [{'H 2 )tp]r is obtained as 



} T P \R 



E II pa* 

{Ji, m } (l,m) 



Z^{Si}\l^(i ) j) J i,j s i^]) e 



+iV B (iV i? - 1) J 2 



ln(r/it) 
ln(p/g) 

P- 9 + 



^{S t } e 

2 

(r + u) 

ln(r/u) ' 
(r — u) 



Mp/q) 

By applying the Cauchy-Schwarz inequality, we obtain 

ln(r/u) 



1(H) 2 



T P \R 



>m) 



T P \R 



N 2 B J 2 



1 2 



p-q + 



\n(p/q) 



(r — u) 



(17) 



(18) 



By using Eqs. (fTUj) . (fT7|) and (|18p. we obtain the rigorous upper bound of the specific 
heat at T = Tp as 



C < k B N B \^ln. 
-2k B N B (ln 



[p + q-{p- qf] + k B N B In 



[r + u 



In 



- q){r - u) . 



(19) 



Eq. (|19|) shows that the specific heat has no singularity when T = Tp. 

It is straightforward to apply the same arguments as in the conventional ± J Ising 
spin glass model case to derive identities and inequalities for correlation functions 
and order parameters. The results are rigorous and are obtained as 



[(SiSjftU = [(SiSjTriSiSjWR (n= 1,3,5,...) 
[P(m)} R = [P(q)] R (T = T P ) , 

IKSiS^rU < [KSiS^TpWn, 
[sgn(S , jS , J ) T ] J R < [sgn(S i S j ) T p}R ■ 



(20) 

(21) 
(22) 
(23) 



We omit the description of deriving Eqs. (f2"0"j) - (|2"3"|) . since these applications are 
straightforward. See Ref. .3). When T = Tp, Eq. (|20p with n = 1 shows that the 
ferromagnetic correlation function on the left-hand side is equal to the spin glass 
correlation function on the right-hand side, and the limit \i — j\ — > oo yields m = q, 
where m and q are the magnetization and the spin glass order parameter respectively. 
The distribution function of the magnetization, [P(m)] R , is 



[P{m)\ 



R 



En* 

{Ji,j} (hi) 



E {Si} S(Nm-EiSi)e^> J ^ 



E 



{Si} • 



(24) 
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and the distribution function of the spin glass order parameter, [P(g)]#, is 

X 1 E^jE^/^^^^ 



(25) 



where the spin Sf J is the spin of the fc-th replica at a site i. Eq. ()2ip shows that the 
distribution of the spin glass order parameter is consistent with the distribution of 
the magnetization when T = Tp. Eq. (|22p shows that the phase boundary between 
the ferromagnetic and non-ferromagnetic phases below the multicritical point should 
be either vertical or reentrant in the phase diagram. Eq. f|23|) shows the spin pair 
becomes mutually parallel ignoring the magnitude when T = Tp. 

The present results described above are given on condition that Eqs. ([3]), (fT2|) 
and (fTUj) are satisfied. The most part of the present results are not the results on the 
so-called Nishimori line, although a part of the present results are equivalent to the 
results on the Nishimori line. The present results are on a three-dimensional solid 
in a four-dimensional space consists of p, q, r and u. If p and q are fixed, Tp/ J is 
also fixed. In addition, even if Tp/ J is fixed, the results are affected by the value 
of a. The value of a depends on the values of r and u even if p and q are fixed. 
If we set p = x/2, q = (1 — x)/2, r = x/2, u = (1 — x)/2 and 1 > x > 1/2, the 
results are results for a = 1 and 2(3pJ = ln[x/(l — x)], and are equivalent to the 
results for the conventional ±J Ising spin glass model, where x is the probability 
that the interaction is ferromagnetic, and 1 — x is the probability that the interaction 
is antiferromagnetic. In addition, if we set p = x/2, q = (1 — x)/2, r = (1 — x)/2, 
u = x/2 and 1 > x > 1/2, the results are results for a = — 1 and 2j3pJ = ln[x/(l— x)], 
and are also equivalent to the results for the conventional ± J Ising spin glass model. 
These results for the conventional ±J Ising spin glass model are results on the 
Nishimori line, and are equivalent to the results in Refs [5]) ,[6]) . If we set r = y/2, 
u = y/2 and 1 > y > 0, the results are results for a = and 2[3pJ = \n(p/q), and are 
equivalent to the results for the bond-diluted ± J Ising spin glass model, where y is 
the probability that the interaction is diluted. The results for the bond-diluted ± J 
Ising spin glass model are equivalent to the results in Ref[7]). Therefore, the present 
solutions are general solutions which include the solutions^® for the conventional 
dbJ Ising spin glass model and the solutions^ for the bond-diluted ±J Ising spin 
glass model. 

An extended ± J Ising spin glass model were investigated by using a gauge sym- 
metry. This model has iJi interactions and ±J2 interactions. We showed that a 
gauge symmetry is usable to study this model. The exact internal energy, the rig- 
orous upper bound of the specific heat and some rigorous relations for correlation 
functions and order parameters were shown by using the gauge symmetry. The re- 
sults are rigorous, and do not depend on any lattice shape. A part of our results, e.g., 
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the value of the exact internal energy should be useful for checking the computer 
programs for investigating this model. In addition, we found that the present solu- 
tions are general solutions which include the solutions® on the Nishimori line for 
the conventional ± J Ising spin glass model and the solutions® for the bond-diluted 
± J Ising spin glass model. 

For this article, it was pointed out that, in Ref. 6), it is already mentioned 
that this distribution (Eq. ([2]) in this article) is applicable for the local gauge 
transformation: this distribution is exactly the same form as Eq. (8) in Ref. [6]) 
(P(Jij) = f(Jij) exp(/3 p Jjj) and f(—Jij) = f(Jij)) which are mentioned as an ap- 
plicable distribution of the local gauge transformation. We think this suggestion is 
correct. On the other hand, in Ref. [S}, the gauge transformation is not explicitly ap- 
plied to this distribution (this model), and, in Ref. [6]), the exact internal energy and 
so forth for this model are not explicitly derived. Probably, this article is the first ar- 
ticle for study of explicit application of the gauge transformation to this distribution 
(this model). 
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